Programmieren 
mit 
FORTRAN 


aus: ComputerKurs-Magazin, Jahrgänge 1984/1985 


Frühe Formeln 


Obwohl FORTRAN schon Mitte der fünfziger Jahre entwickelt wurde, 
hat sich die Sprache viele ihrer ursprünglichen Merkmale erhalten. 
Im ersten Artikel dieser neuen Serie sehen wir uns den allgemeinen 
Programmaufbau an und untersuchen, wie Ein- und Ausgabe und die 
mathematischen Abläufe funktionieren. 


ORTRAN wird ım allgemeinen als die erste 

Hochsprache angesehen. Die Sprache wur- 
de ın einer Zeit entwickelt, als Computer fast 
ausschließlich für die Zahlenverarbeitung ein- 
gesetzt wurden und kaum Möglichkeiten der 
Stnngverarbeitung hatten. FORTRAN ıst noch 
Immer eine der wenigen Sprachen, die kom- 
plexe Zahlen dırekt verarbeiten. Mit ihrer ma- 
thematischen Struktur lassen sıch schnelle und 
leistungsfähige Routinen schreiben, die sıch 
auch für Grafik eignen. Es gibt eine Reihe Ma- 
schinen, deren Betriebssysteme zumindest teil- 
weıse in FORTRAN geschrieben wurden. 

Über Jahrzehnte entstand eine große Zahl 
Subroutinen für ein breites Spektrum numeri- 
scher Änwendungen. Dies ist einer der Haupt- 
gründe für die Beliebtheit derinzwischen schon 
recht betagten Sprache. Da FORTRAN die eiın- 
fachste Methode ist, diese Software einzuset- 
zen, scheint sie noch für einige Zeit eines der 
wichtigsten Programmierinstrumente für Wiıs- 
senschaftler und Ingenieure zu bleiben. 

Auch für Microcomputer gibt es FORTRÄN- 
Versionen, die jedoch nur selten den vollen 
Standard des aktuellen FORTRAN 77 erreichen. 
Die meisten sind Abwandlungen von FORTRAN 
IV mit entsprechenden Erweiterungen. Wır 
werden mit FORTRAN IV beginnen, da diese 
Version den besten Einblick in die Sprache bie- 
tet und später die Erweiterungen und Äbwei- 
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chungen vom Standard erläutern. 

Der Aufbau von FORTRAN wird wesentlich 
von der Tatsache geprägt, daß früher für die 
meisten Ein- und Ausgaben 80-spaltige Loch- 
karten verwandt wurden. Die letzten acht Spal- 
ten einer Karte enthielten normalerweise die 
Kartennummer, so daß für jeden Befehl 72 Spal- 
ten zur Verfügung standen. FORTRAN kann da- 
her pro Zeile nur einen Befehl mit einer Länge 
von maximal 72 Zeichen verstehen. Wenn ein 
Befehl sıch über mehrere Zeilen erstrecken 
soll, müssen die folgenden Zeilen als Fortset- 
zung markiert sein. Die 72 Spalten werden noch 
weiter eingeschränkt: Die ersten fünf Spalten 
sind für die Nummer des Befehls resemnert, 
während die sechste anzeigt, ob eine Fortset- 
zungszeile vorliegt. Der eigentliche Befehl be- 
ginnt erst ab Spalte 7. 

FORTRAN-Zeilen brauchen nicht wie in BA- 
SIC numenert zu werden. Jede Zeile, die von 
einem anderen Befehl angesprochen wird, muß 
Jedoch mit einer Zahl als Label versehen seın. 
Diese Nummern brauchen nicht ın einer be- 
stimmten Reihenfolge stehen, müssen aber eın- 
deutig sein. Ein C ın Spalte Eins markiert die 
Zeile als Kommentar. Weil Datentypen und 
Steuerstrukturen fehlen, sind FORTRAN-Pro- 
gramme schwer zu strukturieren. Da die Spra- 
che schwer zu lesen ist, spielen Kommentare 
eine wichtige Rolle. 
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Wie ın BASIC müssen Varlablen nicht dekla- 
riert werden. Sie lassen sich an jedem Punkt 
des Programms einfach durch Angabe eines 
neuen Namens einführen. Damit tritt die vom 
System vorgegebene Vereinbarung in Kraft, 
daß Varıablennamen, die mit IJ,K,L,M oder N 
anfangen, Ganzzahlen bezeichnen und alle an- 
deren Namen reelle Zahlen. 

Die Komplexität von Ein- und Ausgabe berei- 
tet FORTRAN-Anfängern die meisten Schwie- 


Vergleichende Logik 


Logische Operatoren und Verbindungselemente: 


FORTRAN BASIC 


‚NE. a, 
x6]E e- 
.GE. >= 
Kl. = 


<= 


NOT 
AND | 
OR | 


rngkeiten. Auch dieses Problem entstand aus 
der ursprünglichen Verwendung des starren 
Lochkartenformats. Jeder E/A-Vorgang braucht 
zwei Befehle: READ oder WRITE bestimmen 
das E/A-Gerät und die eingesetzten Varlablen, 
während FORMAT Typ und Anordnung der Da- 
ten auf der Lochkarte, ım Datensatz oder der 
Eingabezeile des Terminals angibt: 
READ(1,100) X, I, ICHAR 

100 FORMAT(F7.2,14,A1) 

READ veranlaßt die Ausführung eines Eingabe- 
vorgangs. Die erste Zahl hinter der Klammer be- 
stimmt das Eingabegerät. Die Gerätenummern 
werden von der Soft- und Hardwareimplemen- 
tatıon festgelegt. Die Nummern bezeichnen 
Terminals, Drucker (nur Ausgabe), Lochkarten- 
leser oder Stanzer und die Dateien der Diskette 
oder Cassette. 

Die zweite Zahl der Klammer bezieht sıch auf 
die Nummer des entsprechenden FORMAT-Be- 
fehls. (Beachten Sıe, daß FORMAT in unserem 
Beispiel mit der gleichen Zahl versehen 
wurde.) Da FORMAT nicht ausgeführt wird, son- 
dern nur Informationen liefert, kann es an einer 
beliebigen Stelle des Programms stehen. Die 
Befehlsnummer ist der einzige Bezug, und so 
können mehrere READs oder WRITESs das glei- 
che FORMAT einsetzen. Der READ-Klammer 
folgt die Vanablenliste, denen die Eingabe- 
werte zugeordnet werden sollen. 

FORMAT enthält eine Beschreibung jeder 
angegebenen Vanable. Hier die wichtigsten: 
@ F: F7.2 zeigt beispielsweise an, daß die Ein- 
gabe für X eine reelle Zahl ist, die ersten sieben 
Spalten der Zeile, Karte oder des Datensatzes 
belegt und hinter dem Punkt zwei Ziffern hat. 
@ I: 14 bestimmt die Eingabe in I als Ganzzahl, 
die vier Spalten belegt. 


® A: Allegt den Wert der nächsten Spalte als 
Zeichen fest. 

FORTRAN kann Zeichen nur verarbeiten, 
wenn sıe als Ganzzahlvariablen gespeichert 
werden. Die maximale Zahl pro Variable hängt 
von der Länge ab, die das System einer Ganz- 
zahlvarlablen zuweist. Bei Ganzzahlen im 16-Bit- 
Formatkönnen nur Alund A2verwandt werden. 
In einigen Sprachversionen ist es jedoch auch 
möglich, dafür Variablen für reelle Zahlen ein- 
zusetzen. Da sıch Varlablen sowohl numerisch 
als auch alphanumerisch verarbeiten lassen 
und die Speicherung nicht immer dem ASCII- 
Code entspricht, kann hier leicht Verwirrung 
entstehen. So müssen längere Strings z. B. über 
Ganzzahlenarrays bearbeitet werden. 

Die eingegebenen Werte müssen genau den 
Vanablenbeschreibungen entsprechen, da 
sonst Äblauffehler auftreten. Die meisten mo- 
dernen Systeme verfügen jedoch über ein 
„freies Eingabeformat", beispielsweise READ 
(1*JAB,C oderauch READ ABC, das sich spe- 
ziell für die Tastatureingabe eignet, bei der ex- 
akte Abstände zwischen den einzelnen Werten 
nur schwer einzuhalten sınd. Dieses freie Ein- 


Schnelle Zahlen 


FORTRAN IV kennt eine Reihe numerischer 
Datentypen: 

@® INTEGER und REAL sind Ganzzahlen und 
reelle Zahlen 

@ DOUBLE PRECISION sind Fließkommazahlen, 
die über die doppelte Menge an Speicher- 
wörtern verfügen wie normale reelle Zahlen. 
®© COMPLEX steht für komplexe Zahlen, die 
einen reellen und einen imaginären Teil 
besitzen. | 
© LOGICAL entspricht dem PASCAL-Typ BOO- 
LEAN und kann die Werte .TRUE. und .FALSE. 
annehmen (die Punkte sind notwendig). 

@ Variablennamen können bis zu sechs Zei- 
chen lang sein. Sie dürfen nur alphanumeri- 
sche Zeichen enthalten. Das erste Zeichen 
muß ein Buchstabe sein. 

® Variablen lassen sich nur am Änfang des 
Programms mit einem Deklarationsbefehl an- 
legen, etwa: 


INTEGER NUM1,NUM2 
LOGICAL FNSH 


gabeformat arbeitet wıe der BASIC-Befehl IN- 
PUT. Ausgaben (mit dem WRITE-Befehl) funk- 
tıonıeren ähnlıch: 


WRITE(1,200) A,B,C 
200 FORMAT(3F7.2) 


Beachten Sıe, daß Beschreibungen sıch mehr- 
fach einsetzen lassen und nicht für Jede Varia- 
ble neu angelegt werden müssen. 

Zeichenketten lassen sıch mit der speziellen 
H-Beschreibung ausgeben. Sie ist jedoch recht 
umständlich, da darın die exakte Zeichenzahl 
aufgeführt sein muß: 


WRITE(1,300)ANS 
300 FORMAT(18HDIE ANTWORT 
LAUTET,F7.2) 
Beı neueren FORTRAN-Versionen lassen sıch 
auch ın AÄnführungsstrichen eingeschlosse 
Strings einsetzen: 
WRITE(1,300)ANS 
300 FORMAT(,DIE ANTWORT 
LAUTET',F7.2) 
Beı der Ausgabe an einen Drucker setzt FORT- 
RAN das erste gedruckte Zeichen zur Steue- 
rung des Papiers eın. Ein typisches System 
kann beispielsweise das Leerzeichen für nor- 
male Zeichenabstände und die Zeichen 12,3 
etc. für die gewünschte Anzahl Leerzeilen ver- 
wenden. Da diese Eigenheit oft vergessen wird, 
geht leicht das erste Zeichen der Ausgabe ver- 
loren, während der Drucker unvorhersehbare 
Dinge mit dem Papıer anstellt. Hier eın Beispiel 
für die Druckerausgabe: 
WRITE(2,400)1,J,K,L,M,N 
400 FORMAT(1H,2(16,5X)) 
Beachten Sıe das zusätzliche IH, für dıe Papıer- 
steuerung. X erzeugt dıe angegebene Zahl 
Leerzeichen. Eine Beschreibungsfolge kann 
geklammert werden — zum Beispiel 2(16,5X). 
Wır verwenden ın allen Beispielen dıe Gerä- 
tenummer ] für dıe Ein- und Ausgabe per Ter- 
mınal. Der korrekte Wert hängt allerdings von 
der eingesetzten Implementation ab. 


Sterne statt Pfeil 


Das mathematische Format und die Zuordnun- 
gen von FORTRAN sınd BASIC sehr ähnlich, mit 
dem einzigen Unterschied, daß FORTRAN für 
Potenzen statt des Aufwartspfeils (t) doppelte 
sterne (**) einsetzt. In den mathematischen 
Ausdrücken lassen sıch reelle Zahlen, Ganz- 
zahlen und sogar Werte mit doppelter Prazısion 
vermischen. Dabeı ıst Jedoch Vorsicht geboten, 
da leicht eın Teil der Berechnung mit der Ganz- 
zahlenarıthmetik ausgeführt wırd, auch wenn 
das Ergebnis eın reeller Wert seın soll. 

=2 

Je 

A=1/J) 

Hier hat A den Wert O, da die Teilung ın der 
Ganzzahlenarıthmetik erfolgte. Die Funktion 
FLOAT kann dıe Ganzzahlenwerte Jedoch ın 
reelle Zahlen verwandeln: 

A=FLOAT(I)/FLOAT(J) 
lıefert das korrekte Ergebnis. 

Da FORTRAN als mathematische Sprache 
konzipiert wurde, besitzt sıe selbstverständlich 
viele Standardfunktionen. Die Liste ist Jedoch 
zu lang, um sıe hier abdrucken zu können. Sıe 
enthält außer den üblichen tnogonometrischen 
und logarıthmischen Funktionen auch viele, dıe 
nur für Mathematiker und Ingenieure von Be- 
deutung sınd. Die meisten Funktionen haben 
eine Version für reelle Zahlen, eine für Zahlen 
mit doppelter Präzision und, falls nötig, auch 
noch das Ganzzahlenformat. 


Einer der Hauptnachteile von FORTRAN war 
das Fehlen von Steuerstrukturen. FORTRAN 77 
beseitigt dieses Problem Jedoch teilweise. So 
wird der vertraute GOTO-Befehl sehr oft einge- 
setzt. Er kann Jeden mit Nummer versehenen 
ausführbaren Befehl (nıcht aber FORMAT) an- 
sprechen. 

Es gibt zweı Vanationen des IF-Befehls: Das 
logische IF entspricht dem IF ın BASIC. Es hat 
das Format: 

IF (logischer Ausdruck) ausführbarer Befehl 
Das arıthmetische IF wırd oft als das schlechte- 
ste Beispiel einer Steuerstruktur zitiert. Es hat 
folgendes Format: 

IF (arıthmetischer Ausdruck) S1, S2, S3 
wobei Sl, S2 und S3 Befehlsnummern sind. Die 
Steuerung wird an die Befehlsnummer S| über- 
geben, wenn der Wert des arıthmetischen Äus- 
drucks kleiner als Null ıst, an S2, wenn er gleich 
Null ıst und an S3, wenn er darüber liegt. Für dıe 
Programmierung eıner dreifachen Verzweı- 
gung Ist dies sehr praktisch, führt aber zueinem 


noch erheblich komplizierteren 


„opaghetti- 


code" als GOTO. 


Durchschnitt 


NINO 


2N22znDnND9 


DD 
rn OD 


DAS PROGRAMM LIEST ZAHLEN 

UND DRUCKT DEN DURCHSCHNITT AUS 
DER INPUT WIRD DURCH DIE EINGABE 
EINER 

NEGATIVEN ZAHL BEENDET 

ZAEHLER UND SUMME INITIALISIEREN 


SUMME=0 
IZAEHLER=0 


NAECHSTE ZAHL LESEN 
READ(1,10)X 

TESTEN, OB SIE NEGATIV IST 

IF (X) 300,200,200 

POSITIVE ZAHL 
IZAEHLER=IZAEHLER+1 
SUMME=SUMME+X 

GOTO 100 

ENDE DER EINGABE 
AVGE=SUMME/FLOAT(IZAEHLER) 
WRITE (I,20)AVGE 

STOP 

FORMATIEREN 

FORMAT (F9.2) 

FORMAT (20HDER DURCHSCHNITT 


IST,F9.2) 
END 
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Zahlenakrobatik 


Das Array ist zwar die einzige Datenstruktur von FORTRAN, doch 
erhält die Sprache viel Flexibilität durch den komplexen Einsatz von 
Subroutinen. Wir untersuchen die speziellen Techniken für die 
Zahlenverarbeitung mit Hilfe dieser Sprache. 


as Array Ist die einzige Datenstruktur, die 

dem FORTRAN-Programmierer zur Verfü- 
gung steht. Wie ın BASIC müssen Arrays vor 
ihrem Einsatz deklariert werden — direkt am 
Programmanfang, da FORTRAN eıne compi- 
lierte Sprache ist. Diese Deklaration erfolgt mit 
dem Befehl DIMENSION, der die gleiche 
Funktion hat wie DIM in BASIC. 

Bei systemvorgegebenen Datentypen (alle 
Varnablennamen, die mit1J,K,L,M oder N an- 
fangen, stellen Ganzzahlen dar und alle ande- 
ren Namen reelle Zahlen) reserviert 

DIMENSION A(20), I(50) 


den Speicherplatz für ein Array mit bis zu 20 
reellen Zahlen und ein weiteres Array für bis 
zu 50 Ganzzahlen. Befehle, die dıe Systemein- 
stellung verändern (INTEGER, REAL, DOUBLE 
PRECISION, COMPLEX oder LOGICAL) müs- 
sen zusätzlich zu der Anweisung DIMENSION 
gegeben werden. Arrays lassen sich aber auch 
gleich ın der Typendeklaration definieren: 

INTEGER ARR1(20) 

REAL NUMS(30) 


Der erste Befehl reserviert Platz für ein Array 
mit 20 Ganzzahlen und der zweite für 30 reelle 
Zahlen. 

Zweidimensionale Arrays sind möglich, ei- 
nige Sprachversionen erlauben sogar drei 
oder mehr Dimensionen. Zweidimensionale 
Arrays werden auf übliche Weise deklarıert: 

DIMENSION ARR(20,30) 

INTEGER ARR(20,30) 


Jede Ganzzahlkonstante oder Varlable läßt 
sıch als Arrayindex einsetzen, beispielsweise 
ARR(3,4) oder INTARR(I). 

FORTRAN IV verarbeitet zwar keine Zei- 
chenstrings, doch stehen bei FORTRAN 77 da- 
für andere Möglichkeiten zur Verfügung (mit 
FORTRAN 77 werden wir uns ın der nächsten 
Folge genauer beschäftigen). Zeichen werden 
wie Ganzzahlen ın normalen Ganzzahlvaria- 
blen gespeichert. In der letzten Folge haben 
wir gesehen, wıe sich ein oder mehr Zeichen 
in Ganzzahlvaraablen unterbringen lassen. 
Zwar kommt das FORTRAN-Ganzzahlenarray 
einem Zeichenstring sehr nahe, doch ist es im- 
mer noch recht unbefriedigend, wenn man die 
Stringverarbeitungsmöglichkeiten von BASIC 
gewöhnt ist. 


Außer den Befehlen IF und GOTO ist die 
DO-Schleife die wichtigste Steuerstruktur von 
FORTRAN. Sie funktioniert wie die FOR... 
NEXT-Schleife in BASIC und hat folgendes 
Format: 

DO Sn (intvar) = Startwert, Endwert, 

Schrittwert 

Sn (letzter Befehl des zu wıederholenden 

Blocks) 


Sn Ist die Nummer des letzten Schleifenbe- 
fehls. Schleifenzähler (intvar) kann jede Ganz- 
zahlvarıable sein. Die Anfangs-, End- und 
Schrittwerte sınd angegeben. In BASIC würde 
die DO-Schleife so aussehen: 
FOR ıntvar = Startwert TO Endwert STEP 
Schrittwert 


NEXT intvar 
Folgende Befehle setzen eın Array auf Null: 
DO 100 I = 1,20, 1 


100 ARR(I)=0 
Der Schnittwert Eins darf weggelassen werden: 
DO 100 | = 1,20 
100 ARR(I)=O 


Die numenierte Anweisung, die das Schleifen- 
ende kennzeichnet, kann jeder ausführbare 
Befehl außer IF, GOTO oder DO sein. 
Manchmal ist es praktisch, keine feste 
Schleifenendanweisung zu haben, da sich der 
Code damit komprimieren läßt. Diese Möglich- 
keit kann aber auch zu schweren Fehlern füh- 
ren, und besonders bei verschachtelten 
Schleifen werden Programme dann völlig un- 
verständlich. Die folgende Initialisierung eines 
zweidimensionalen Arrays Ist jedoch legal: 


DO 1001 = 1,20 
DO 100.) = 1,80 
100 ARR(I,J)=0 


(Beachten Sıe, daß bei beiden Schleifen die 
gleiche Endanweisung verwendet wird.) 

Oft wırd der „Leerbefehl“ CONTINUE einge- 
setzt, um das Schleifenende zu markieren und 
den Code einfacher zu gestalten. CONTINUE 
ist ein FORTRAN-Befehl, der keine Wirkung 
hat und an jeder beliebigen Position eingesetzt 
werden kann. Jede DO-Schleife muß normaler- 
weise mit einem eigenen CONTINUE abge- 


schlossen seın. Das vorige Beispiel sieht dann 
in ausführlicher Schreibweise so aus: 
DET I 12 
DO 100 J = 1,30 
ARR(I,J)=0 
100 CONTINUE 
101 CONTINUE 


Damit FORTRAN-Programme generell über- 
sichtlicher werden, verlangen einige Program- 
mierer sogar, daß alle DOs, IFs und GOTOs mit 
einem CONTINUE abgeschlossen werden. 
Das speziell für Mathematik ausgelegte 
FORTRAN besitzt Anweisungen, die beson- 
ders gut für den Umgang mit Zahlenarrays ge- 
eignet sind. So kann mit einem einzigen 
READ- oder WRITE-Befehl ein ganzes ein- 


oder zweidimensionales Array gelesen oder 
geschrieben werden, wobei FORMAT den Auf- 
bau jeder Zeile oder jedes Datensatzes angibt. 
Zum Beispiel: 

DIMENSION IARR1(20) 


WRITE(1,10)IARRI 


veranlaßt die Ausgabe des gesamten Inhalts 
von IARRI (20 Ganzzahlen). 

10 FORMAT (2014) 
erzeugt 20 vierstellige Ganzzahlen in einer 
Zeile; 

10 FORMAT (14) 


erzeugt eine Ganzzahl per Zeile; 
10 FORMAT (514) 


Wurzelziehen 


EIN FORTRAN PROGRAMM, DAS SETS VON 
DREI WERTEN A,B;C, LIEST, DIE DIE 
KOEFFIZIENTEN DER QUADRATISCHEN 
GLEICHUNG A*X**2+B*X+C=0 ÜBER DIE 
ÜBLICHE FORMEL DARSTELLEN 
FUNCTION DISCR(A,B,C) 
DISCR=B*B-4.0*A*C 
RETURN 
END 
SUBROUTINE SOLVE(A,B,C,X1,X2,NROOT) 
D=DISCR(A,B,C) 
ZWAR LIESSE SICH HIER EIN 
ARITHMETISCHES IF 
EINSETZEN, DOCH IST ES NICHT 
EMPFEHLENSWERT 
IF (D.GE.0.0)GOTO 100 
C DA D UNTER NULL LIEGT, IST DIE WURZEL 
C KOMPLEX 
NROOT=0 
X1=0.0-B/(2.0*A) 
X2=SOR(ABS(D)) 
GOTO 300 
100 IF (D.GT.0.0)GOTO 200 
C _DIST NULL, DAHER GLEICHE WURZELN 
NROOT=1 
X1=0.0-B/(2.0*A) 
KON 
GOTO 300 
C _DIST POSITIV, DAHER ZWEI WURZELN 
NROOT=2 
200 X1=(0.0-B+SORT(D))/2.0*A 
X2=(0.0-B-SORT(D))/2.0*A 
300 RETURN 
END 
C ERST DIE ANZAHL DER SETS LESEN 
C _READ(2,10)NVAL 
C _UÜBERSCHRIFTEN ANZEIGEN 
WRITE(3,11) 
DO 500 L=1, NVAL 
C _DREIWERTE LESEN UND PRÜFEN, OB SIE 
C GEEIGNET SIND 
100 READ (1,12) A,B,C 
IF (A.EQ.0.0) GOTO 500 
C DIE WERTE SIND GEEIGNET. DAHER 


euslsZs 


DIL 


C 


SUBROUTINE AUFRUFEN 

CALL SOLVE(A,B,C,X1,X2, NROOT) 
IF (NROOT.GT.0)GOTO 200 
KOMPLEXE WURZELN 

WRITE (1,13) A,B,X1,X2 

GOTO 500 


200 IF (NROOT.GT.1)GOTO 300 


& 


C 


GLEICHE WURZELN 
WRITE (1,14)A,B,C,X1,X2 
GOTO 500 

ZWEI WURZELN 


300 WRITE (1,15)A,B,C,X1,X2 


» 


GOTO 500 
UNGÜLTIGE WERTE FÜR A, BUND C 


400 WRITE (1,16)A,B,C 
500 CONTINUE 


G 
10 
11 


15 


16 


STOP 

FORMAT BEFEHLE 

FORMAT (14) 

FORMAT(1H ‚8X,1HA,8X,1HB,8X,1HC,8X, 
AHTYPE,8X,2HX1,8X,2HX2) 

FORMAT(1H ‚3F7,2) 

FORMAT(1H ‚4X,3(F7.2,2X),4X, 7HCOMPLEX, 
X,F7.2,3X,F7.2) 

FORMAT(1H ‚4X,3(F7.2,2X),4X,5HEQUAL3X, 
F7.2,3XF7.2) 

FORMAT(1H ‚4%X,3,(F.2.2X).4X, 
7HUNEQUAL3,X, 

F7.2.3X,F7.2) 

FORMAT(1H ‚4X,3,(F7.2,2X),4X, 7HINVALID) 
END 


Durchschnittsberechnung 


199209 


EINE FORTRAN-FUNKTION, DIE 

DEN DURCHSCHNITT 

(IN REELLEN ZAHLEN) EINES STETS MIT N 
GANZZAHLEN 

DES GANZZAHLENARRAYS IARR BERECHNET 
FUNCTION AVGE (lIARR,N) 

DIMENSION IARR (N) 

SUM=0 

DO 100 I=1,N 

SUM=SUM+FLOAT (IARR(I)) 


100 CONTINUE 


AVGE=SUM/FLOAT (N) 
RETURN 
END 
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Da FORTRAN keinen 
fest definierten Befehl 
für das Schleifenende 
besitzt, können mehrere 
Schleifen mit der glei- 
chen Endanweisung ab- 
geschlossen werden. 
Damit läßt sich sehr 
kompakter Programm- 
code schreiben (siehe 
Bild). In BASIC ist der 
gleiche Ablauf zwar 
weniger komprimiert, 
aber dafür auch weit 
leichter zu lesen. 


100 FOR I=1 TO 20 
110 FOR J=1 TO 30 
120 LET A(1,J)=0 
130 NEXT J 

140 NEXT I 


Das letzte 
Element 


vier Zeilen mit je fünf Zahlen und 
10 FORMAT (20A1) 


einen String von 20 Zeichen in eıner Zeile. 
Beim Einlesen oder Ausgeben von zweiıdi- 
mensiıonalen Arrays entsteht oft ein Problem, 
da FORTRAN, ım Gegensatz zu anderen Spra- 
chen, zuerst die Spalten und dann erst die Zeı- 
len speichert. Wenn Sıe hier nicht aufpassen, 
entsteht eine versetzte Version des Arrays. Es 
gıbt Jedoch eine spezielle Form der DO- 
Schleife, die „implizierte Schleifenanweisung"“, 
die sich gut für die Arbeit mit einem Teilarray 
eignet. Implizierte Schleifenanweisungen wer- 
den mit dem Variablennamen ın die Befehle 
READ oder WRITE eingefügt. Die folgenden 
Zeilen füllen beispielsweise das Stringarray 
IARR2( 10,20) mit zehn Strings zu Je 20 Zeichen: 
DO 100 I=1,10 
READ (1,10) (IARR2(I,J), J=1,20) 
10 FORMAT (20A1) 
100 CONTINUE 


FORTRAN bietet die Möglichkeit, mıt Subrouti- 
nen zu arbeiten — wenn auch mit einigen Eın- 
schränkungen. Die Subroutinen sınd völlig un- 
abhängig und arbeiten hauptsächlich mit loka- 
len Vanablen. Sie werden entweder an das 
Ende des Hauptprogramms angehängt oder 
separat geschrieben und compiliert. Subrouti- 
nen können durchaus die gleichen Zeilennum- 
mern einsetzen, dıe auch das Hauptprogramm 
verwendet (in den meisten Fällen sogar die 
gleichen Vanablennamen). Einfache Parame- 
ter werden durch Wertaufruf übergeben, wo- 
beı Werte, dıe von der Subroutine geändert 
wurden, auch die entsprechenden Variablen 
des aufrufenden Programms beeinflussen. 
Globale Varlablen, auf dıe sıch Subprogramme 
und Hauptprogramm beziehen, müssen ın beı- 
den Programmodulen mit COMMON dekla- 
niert werden. Es lassen sıch eine ganze Reıhe 


i 100 ARR N 


(1, = 0 
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von COMMON-Blöcken anlegen und mit Na- 
men versehen, doch meistens genügt eın ein- 
zıger Block dafür. 

Die Vanablen der COMMON-Anweisungen 
müssen einander ın Typ und Gesamtlänge ent- 
sprechen. So könnte das Hauptprogramm beı- 
spielsweise den Befehl 

COMMON 11,12,13 


enthalten und dıe Subroutine 
COMMON |I(3) 


In beiden Fällen wird der Speicherplatz auf 
die gleiche Weise eingesetzt, wobeı dıe Para- 
meter sıch ım ersten Fall auf Ganzzahlenvaria- 
blen beziehen und im zweiten Fall auf drei Ar- 
rayelemente des gleichen Typs. 

Ein typischer Subroutinenaufruf sieht folgen- 
dermaßen aus: 

CALL SUBNAME(X,Y,l,J) 


wobei die Subroutine mit dem Befehl 
SUBROUTINE SUBNAME(A,B,L,M) 


anfängt. (Beachten Sıe, daß dıe Parameter ım 
Inneren der Klammern ın Bezug auf Anzahl 
und Typ übereinstimmen.) 

Die Anweisung RETURN gibt die Steuerung 
an das aufrufende Programm zurück (in einer 
Subroutine können sıch mehrere RETURNSs be- 
finden). Eines der RETURNs sollte Jedoch in 
jedem Fall ausgeführt werden — unabhängig 
von dem Weg, der durch die Subroutine führt. 
In Subroutinen wie ım Hauptprogramm muß 
der letzte Befehl ımmer eın END seın. 

Bei der Parameterübergabe ist es in be- 
stimmter Hinsicht möglich, die Größe eınes Ar- 
rays dynamisch zu verändern. Obwohl das ın 
einer Interpretierten Sprache wie BASIC leicht 
auszuführen ist, gibt es diese Möglichkeit ın 
compilierten Sprachen nur selten, da der Speiı- 
cherplatz bei der Compilierung fest zugeord- 
net wırd und nicht erst mit Geschwindigkeits- 
verlust beim Ablauf des Programms. 

Da Subroutinen nur einen einzigen Wert zu- 
rückliefern, können sıe auch (wıe ın PASCAL) 
als FUNCTION angelegt werden. Derartige 
Subprogramme werden, wie normale Subrouti- 
nen, separat geschrieben, beginnen aber mit 
der Anweisung FUNCTION. Vor dem RETURN 
sollte eine Zuweisung auf den Funktionsna- 
men liegen. Der Funktionsname selbst be- 
stimmt den von der Funktion zurückgelieferten 
Datentyp, wobei dıe Vorgaben der vom System 
vorgegebenen Varlablennamen angewandt 
werden. Diese Anforderung läßt sıch zwar um- 
gehen, sollte aber wenn möglıch beibehalten 
werden. Das Hauptprogramm setzt die defı- 
nıerten Funktionen auf die gleiche Weise eın 
wıe die normalen Systemfunktionen. Der Funk- 
ttıonsname (mit zugehöriger Parameterliste) 
kann dabeı selbstverständlich ın jedem Aus- 
druck erscheinen, in dem eıne Vanable dieses 
Typs auftreten darf. 


In den meisten BASIC- 
Versionen sind ver- 
schachtelte IFs mog- 
lich. FORTRAN hat je- 
doch zusätzlich die 
Struktur ELSEIF, mit der 
sich IF-Anweisungen 
uber mehrere Pro- 
grammzeilen bis zu 
einer beliebigen Tiefe 
verschachteln lassen. 
ENDIF beendet die 
Struktur. 
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Aus alt mach neu 


Da FORTRAN aus den fünfziger Jahren stammt, fehlen dieser Sprache 
viele Strukturen und Möglichkeiten, die in modernen 
Programmiersprachen als selbstverständlich angesehen werden. Wir 
sehen uns einige Techniken an, die diesen Mangel beheben. 


Eine sehr praktische Form dıeser Änweı- 
sung arbeitet mıt verschachtelten IF's: 


ORTRAN fehlen zweı wıchtige Eıgenschaf- 
ten: Es hat keine brauchbaren Steuerstruk- 
turen, die die Sprache leıcht verständlich ma- 
chen und verfügt außerdem nur über eınge- 
schrankte Möglıchkeiten der Zeıichenverarbeı- 
tung. Im Laufe der Zeit gab es eıne Reıhe von 
Versuchen, diese Einschränkungen ganzlıch 
zu beseitigen. 
Eine Möglıchkeit dafür bietet der Eınsatz 
eines „Vorprozessors" (eine Art Compiler auf 
einer hoheren Ebene). Er bearbeıtet Pro- 


IF (logischer Ausdruck) IHEN 


gramme, dıe mit einer erweiterten FORTRAN- ELSEIF 

Version erstellt wurden und setzt für alle Zu- 
satzeigenschaften dıe entsprechenden Ab- Else 

laufe des Standard-FORTRAN eın. Damit las-- 0.0... 

sen sıch struktunerte Programme schreiben, ENDIF 


Es lassen sıch beliebig viele ELSEIFSs für tıe- 
fere Strukturierung verwenden. 


dıe ın kompilierter Form dem Standard ent- 
sprechen und sıch allgemein verwenden las- 
sen. WATFOR und RATFOR sınd dıe am weiıtest 
verbreiteten Versionen. 

Die Eigenschaften von FORTRAN 77 ent- 
sprechen denen der bekannteren Vorprozes- 
soren. Aus diesem Grund blıeb dıe Essenz der 
Sprache erhalten, obwohl vıele Compiler für 
Mıcrocomputer auf dem erweiterten FORTRAN 
IV aufbauen und nicht auf FORTRAN 77. 

FORTRAN besıtzt keine Blockstruktur wie 
ALGOL, PASCAL oder C, was sıch speziell beı 
der Einrichtung von Steuerstrukturen negativ 
bemerkbar macht. Außerdem können nıcht 
mehrere Befehle zu einer Anweisung zusam- 
mengefaßt werden (beispielsweise durch BE- 
GIN...END). Blöcke lassen sıch nur ısolieren, 
wenn sie als separate Subroutinen angelegt 
werden, oder von GOTO-Befehlen umgeben 
sınd. Es gıbt Jedoch eıne neue Steuerstruktur, 
dıe allerdings nicht so recht zu der Sprache 
Baht: F,.-TIIEN& ELSE.» ENDIF 


FALSCH 


IF (logıscher Ausdruck) THEN 


FALSCH 


Der logısche Ausdruck hat das gleiche Format 
wıe Standard-FORTRAN. Normalerweise wırd 
der ELSE-Teıl weggelassen. 


Die zweite wesentliche Verbesserung ıst dıe 
Einführung eınes Datentyps für Zeichen. Zeı- 
chenstrings werden — wıe dıe anderen Dekla- 
ratıonen — am Programmanfang ın eınem der 
folgenden Formate festgelegt: 


Dıe maxımale Länge eınes Strings wird hıer 
mit *n angegeben, wobeı n eine Ganzzahl ıst. 
Dieser Vorgang kann entweder auf das Schlüs- 
selwortt CHARACTER angewandt werden 
(wenn alle Strings die gleiche Länge haben), 
oder auf einzelne Strings mit individuellen Be- 
zeichnungen. Die beiden aufgeführten Dekla- 


CODEKNACKER 

C EIN PROGRAMM, DAS GEHEIMNACH- 

C RICHTEN DECODIERT. BEACHTEN SIE DIE 
C _ ANWEISUNG ‚PROGRAM' IN FORTRAN 77 
C 


INTEGER COUNT,PTRIN,PTROUT 
LOGICAL FILE 

CHARACTER *10 NUM, CHAR 
CHARACTER *30 IN, OUT 

DATA COUNT /0/ 


‚DATA ANWEISUNG ZUR INITIALISIERUNG 
EINER VARIABLE 


WENN KEINE GEHEIMDATEI EXISTIERT, 
FEHLERMELDUNG AUSGEBEN 


DSL DF TS 


INQUIRE (FILE="GEHEIM’,EXIST=FILE) 
IF (.NOT.FILE) THEN 

WRITE(1,10) 

STOP 
ENDIF 


DATEIEN ERÖFFNEN 


DEN 


OPEN (UNIT—12,FILE="MELDUNG’,STATUS 
=-'NEW') 


DECODIERSCHLUSSEL ANGEBEN 


CILCHE) 


NUM-'0123456789' 
CHAR ="MW3 OD$% T' 


GEHEIME MELDUNG LESEN UND 
DECODIEREN 


snrnn2p9 


0 READ (11,20, END=1000)IN 
OUT=- ° 
PTROUT-1 
DO 2001--1,30 
PTRIN=INDEX (CHAR.IN(I:I)) 


STANDARDFUNKTION INDEX’ 
BESTIMMT DIE POSITION 
DES TEILSTRING IN(I:1) IM STRING CHAR 
IF (PTRIN.NE.O) THEN 
OUT (PTROUT:PTROUT)=NUM 
(PTRIN:PTRIN) 
PTROUT=PIIROL TER 1 
ENDIF 


DEE 


ratıonen ergeben zweı Strings der Länge 9, 
eınen mit Länge 7 und eınen mit Lange 16. 
Stringarrays werden wıe üblıch definiert: 


deklanert eın Array von 50 Strings mıt Jeweils 
vıer Zeichen. Stnngkonstanten werden String- 
varıablen folgendermaßen zugeordnet: 


Ist der zugeordnete Strnng zu kurz, werden 
Leerzeichen angehängt; ıst er zu lang, werden 
die Zeichen am Stringende ıgnornert. 

Strings können mit normalen relatıonalen 


200 CONTINUE 
WRITE (12,20) OUT 


BEACHTEN SIE, DASS EIN- UND AUSGABE 
DAS GLEICHE FORMAT HABEN 


eo 


COUNT=COUNT+1 
GOTO 100 


C 
C _ENDVORGANG 

C MIT 'ENDFILE’ WIRD DIE MARKIERUNG 
C DES DATEIENDES GESCHRIEBEN 

€ 


1000 ENDFILE (UNIT = 12) 
STOP 


FORMATS 

FORMAT (1H,’ES EXISTIERT 
KEINE GEHEIMDATEI’) 
FORMAT (A) 


sn 


oO 


BEACHTEN SIE, DASS NICHT 
ANGEGEBEN WERDEN MUSS, 
WIEVIEL ZEICHEN 
EINZUGEBEN SIND 


ADITIDTIIO-D 


END 


IT 
Kup 


Dieses FORTRAN-Pro- 
gramm liest Daten aus 
einer Datei namens GE- 
HEIM und dechiffriert 
über einen Decodier- 
schlüssel den chiffrier- 
ten Text. 


FT 
Hi 
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Da FORTRAN-Compiler 
recht teuer sind, gibt es 
sie normalerweise nur 
für kommerzielle Ma- 
schinen. Das im Bild 
gezeigte Paket läuft auf 
dem Tandy 2000. Durch 
die begrenzten Spei- 
cherkapazitäten sind 
FORTRAN-Anwendun- 
gen auf Micros proble- 
matisch, wenn nicht 
unmöglich. 

Mit der wachsenden 
Verfügbarkeit von 
CP/M auf Heimcompu- 
tern wie dem Commo- 
dore 128, Schneider 
CPC und anderen, ste- 
hen vielen Anwendern 
nun auch FORTRAN- 
Versionen zur Verfü- 
gung. Für Besitzer von 
Schneider Geräten gibt 
es eine besonders 
preiswerte Version, das 
„Nevada FORTRAN“. 
Sie kann bei New Star 
Software Ltd., 45 Plo- 
vers Mead, Wyatts 
Green, Essex, CM15 
OPS bestellt werden. 
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Operatoren wıe GT.,.LT etc. verglichen werden. 
Es gıbt Jedoch zwei neue Operatoren nur für 
Stringfunktionen. Der Teilstringoperator lıefert 
eine Zeichenfolge des Strings. 


CHN8:8) 


holt den Teilstrnng von CHI], der beim dritten 
Zeichen beginnt und beim fünften endet. 


CHI="ABCDEFGHI' 
CH2=CHi3’5) 


legt ın CH2 den Inhalt ‚CDE’, oder genauer 
CDE "ab. Beachten Sie, daß die zweite Zahl 
die Position des letzten Zeichens angıbt und 
nicht die Länge des Teilstrings. 

Der Verkettungs-Operator (//) verbindet 
zwei Strings miteinander. 

CHARACTER CH1*4,CH2*4,CH3*8 

CH1='ABCD’ 

CH2="WXYZ 

CH3-CHNrGH2 


ergibt CH3 mit dem Inhalt ABCDWXYZ'. 

Um FORTRAN nnchtig nutzen zu können, Ist 
eın System mit Diskettenlaufwerken nötig. Die 
meisten FORTRAN-Versionen können sequen- 
tiell und direkt auf Dateien zugreifen. Die 
Quelle oder die Bestimmung eınes E/A-Vor- 
gangs wird von der Gerätebezeichnung des 
READ- oder WRITE-Befehls bestimmt. Einige 
Parameterwerte — normalerweise die niedri- 
geren — sind für E/A-Geräte wie Tastatur, Bild- 
schirm und Drucker reserviert, die anderen 
Nummern stehen dem Programmıerer zur Ver- 
fügung. Der Befehl OPEN ordnet eıner Disket- 
tendateı eine Namen zu: 


OPEN (UNIT=Ganzzahlausdruck, 
FiILE=Dateiname, STATUS=Status) 


Der Ganzzahlenwert für UNIT bezeichnet ın 
den folgenden READ- oder WRITE-Befehlen 
die Dateinummer. Status kann mehrere Werte 


annehmen -— beispielsweise OLD für eine Ein- 
gabedatei, die es bereits gıbt, oder NEW für 
eine Ausgabedateı, dıe angelegt werden soll. 
Falls die angesprochene Dateı vom einfa- 
chen sequentiellen Format abweicht, kann 
OPEN um eine Reihe von Zusatzangaben er- 
gänzt werden: ACCESS= bestimmt den direk- 
ten oder sequentiellen Zugriff; FORM= gibt an, 
ob die Datei formatiert oder unformatiert ist; 
IOSTAT bietet die Möglichkeit, Fehler abzufan- 
gen, falls keine Dateizuordnung zustande 
kommt; RECL= legt dıe Dateilänge fest. 
CLOSE hat einen ähnlichen Aufbau: 


CLOSE (UNIT=Ganzzahlausdruck, 
STATUS=Status) 


Der Befehl schließt eröffnete Dateien, wird 
aber nur selten benötigt, da bei Beendigung 
des Programms alle Dateien automatisch ge- 
schlossen werden. 

Von den anderen Dateibefehlen ıst INQUIRE 
möglicherweise am ınteressantesten, da sıch 
damit der aktuelle Status Jeder beliebigen Da- 
tei feststellen läßt. 

Auch zu READ und WRITE sınd Zusätze 
möglıch, die die E/A-Vorgänge bestimmen. 


READ(7,20,REC=I)A,B,C 


liest beispielsweise Record I eıner Datei, die 
für den direkten Zugriff eröffnet wurde. 


READ(7,11,END=1000)A,B,C 


übergibt die Steuerung an Anweisung Num- 
mer 1000, wenn ın einer sequentiellen Dateı 
das Dateiende gefunden wird. 

Zum Schluß noch ein Blick auf die FOR- 
TRAN-Compiler. Sıe arbeıten schnell und ef- 
fektiv und erzeugen schnellen, kompakten 
und leistungsfähigen Code, der besonders beı 
Echtzeitanwendungen große Vorteile bringt. 

Fehler lassen sıch mit den FORTRAN-Com- 
pilern Jedoch nicht gut abfangen. So akzeptie- 
ren die Compiler Fehler, dıe beispielsweise 
vor einem PASCAL-Compiler nie bestehen 
könnten. So haben Leerzeichen in einer FOR- 
TRAN-Anweisung etwa keine Bedeutung. Ein 
Compiler nımmt daher oft als erstes alle Leer- 
zeichen aus einer Befehlszeile. Nun legt der 
Befehl 


DO 100 I=1,100 


eine Schleife an, die 100 mal ausgeführt wırd. 
Eın Punkt anstelle des Kommas ergibt nach Be- 
seitigung aller Leerzeichen Jedoch 


DO100I=1.100 


und damit eine völlig legitime Wertzuweisung 
für dıe Variable DO100I. Ein Gerücht besagt, 
daß dieser kleine Irrtum im Weltraumpro- 
gramm der Vereinigten Staaten zu einem sehr 
teueren Fehler geführt haben soll. FORTRAN 
war lange Zeit die Standard Programmierspra- 
che der Streitkräfte. Fehler dieser Art waren 
Anlaß für die Entwicklung von ADA. 


